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We measure evolution of spectra, spatial correlation functions and probability density functions 
(PDF) of waves appearance for one- dimensional generalized Nonlinear Schrodinger equation: (1) 
accounting for six- and eight- wave interactions terms and (2) accounting for six- wave interactions, 
dumping (linear dissipation and three-photon absorption) and pumping terms. All additional terms 
beyond the classical NLS equation are small. We observe strongly non-Gaussian PDFs with "fat 
tails" in the region of large amplitudes when higher waves appear more frequently. For generalized 
NLS equation with six- wave interactions, dumping and pumping terms we demonstrate absence 
of non-Gaussian addition to PDF for zeroth six-wave interactions coefficient and increase of non- 
Gaussian addition with six- wave interactions term. 



1. Since the first observation by Solli et at. in 2007 [1 , optical rogue waves - large wave events 
with statistics drastically different from Gaussian that appear randomly from initially smooth 
pulses - has drawn much scientific attention from both optical and hydrodynamic society. In case 
of optics rogue waves are huge waves that can damage optical systems and therefore their ap- 
pearance must be controlled. For hydrodynamics optical rogue waves are interesting phenomenon 
that can be conveniently studied in laboratory conditions and that occurs in systems described 
by the similar equations and has similar to hydrodynamic rogue waves statistics [2j . The current 
study of optical rogue waves went in two main directions: 1) harnessing and control of rogue 
waves emergence by seeding the initial pulse [31 and 2) understanding physical mechanisms 
underlying the phenomenon by establishing connection between different linear and nonlinear 
terms in the equations of motion and appearance of non-Gaussian tails in the probability of large 
waves occurrence [5*-^. The aim of the current publication is to study the influence of higher 
nonlinearity, namely six-wave interactions term, on the frequency of rogue waves appearance for 
two different nonlinear systems one of which is conservative and the other is non-conservative. 

Our research can also be considered in the broader context. Nonlinear Schrodinger equation 

i*i + *,, + |*|2* = 0, (1) 

generalizations of which are studied in the framework of rogue waves phenomenon, describes 
many physical systems from Bose-Einstein condensate and plasma oscillations to propagation 
of waves in optics and hydrodynamics. We consider the first nonlinear term beyond Eq. ([T]) 
that naturally appears as next-order term in perturbation theory and investigate its influence 
on the map of extreme events of the system in the regime when dynamics of the system is close 
to that described by the classical NLS equation. In this sense our choice of parameters is more 
natural than in the recent publication [9 where the similar research was made for the quintic 
NLS equation. 

Let us suppose that the current state of a system consists of multitude of uncorrelated linear 
waves, 

^ = ^ a/e exp(A:x - ujkt + ^/c). 

k 

If a/c and (j)k are random uncorrelated values and number of linear waves is large enough, then 
probability density function (PDF) for amplitude |^| is Gaussian-distributed, 



PI)F(I^I) -exp(-|^|V2a2). 
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In the current publication we search for deviations from Gaussian distribution at large amplitudes, 
namely for signs of non-exponential tails when higher waves occur more frequently. Presence of 
such a tail would mean the probability of occurrence of a large wave event is by several orders of 
magnitude higher than that predicted by linear theory and therefore the map of extreme events 
is qualitatively different from the ordinary Gaussian. 

In this work we perform numerical simulations of wave field ^ evolution in the framework of two 
different nonlinear equations for ensemble of 10000 initial distributions ^|^=o for each nonlinear 
system. Inside each ensemble initial distributions differ only by realization of stochastic noise with 
a fixed statistical properties. Based on the ensembles, we measure spectra Ik = (l^/cP) where (..) 
stands for averaging over ensemble, spatial correlation functions g{x) = (^(?/, + x, t)) and 

PDFs for amplitudes of waves |^| and examine how these functions depend on time and different 
nonlinear terms included in the equations of motion. 

The paper is organized as follows. Next section gives brief overview of nonlinear systems we 
examine and also numerical methods we use. Results of our numerical simulations and their 
comparisons are listed in Section 3. Section 4 contains conclusions and acknowledgements. 

2. In this work we consider two different one-dimensional generalized Nonlinear Schrodinger 
equations: 1) accounting for six- and eight- wave interactions terms, 

i'^t + "^xx + 1^1'^ + - = 0, a, /3 > 0, a, p/a < 1, (2) 

and 2) accounting for the collapsing six-wave interactions term, as well as dumping terms (linear 
dissipation and three-photon absorption) and a pumping term, 

z^t + (l-idi)^^^ + + (a + zdn)!^!"^^ = ip^, a,dudn,p > 0, a,di,dn,p<^ 1. (3) 

Here t is time, x is spatial coordinate, a and f3 correspond to six- and eight- wave interactions 
terms respectively, coefficients di and dn determine linear and nonlinear dumping respectively 
and coefficient p corresponds to deterministic forcing of the system. Eq. ([2])-([3| correspond to 
anomalous group velocity dispersion (GVD) regime because optical rogue waves are huge pulses 
lying entirely in the anomalous GVD [T]-[8]. Eq. ([2| is the Hamiltonian one, 

SH 

^^^ = ^' 

with Hamiltonian 

H = Hd + Hi + He + Hs, Hd = J l^.l'dx, 

Hi = - j^-^dx, He = -J dx, Hs = j ^ dx. (4) 

In addition to Hamiltonian, Eq. Q also conserves wave action N = j \^\^ dx and momentum 
P = {i/2) J (^*^ — ^a?^*) dx. ForEq. ^ energy Hq, wave action and momentum 

become functions of time. 

Addition of six-wave interactions term only to the classical NLS equation ([T]) results in gen- 
eration of blow-up collapses in a finite time. There are two ways to regularize these collapses: 
by adding next-order nonlinear term in the form of defocusing eight-wave interactions ([2| or 
with the help of dumping terms Thus, nonlinear dissipation in the form of three-photon 
absorption prevents formation of waves with too high amplitude, while linear dissipation may 
originate from optical filtering JO] and prevents appearance of too high gradients. Energy loss 
due to the dumping terms that manifests itself mainly as a result of regularization of collapses 
(compare to [9 ) is compensated by inclusion of forcing term. In this work we limit ourselves with 
consideration of only deterministic forcing in Eq. ([3| because according to ^ we do not expect 
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qualitative difference for other types of forcing. Specific values of coefficients for Eq. 

a, P^di^dn^p were chosen to be small enough in order to make the dynamics of the systems close 

to that described by the classical NLS equation ([T]). 

We solve Eq. ([2|-(|3| numerically in the box — IGtt < x < IGtt with periodical boundary 
conditions starting from the initial data ^|t=o = 1 + e(x), where e{x) is stochastic Gaussian- 
distributed noise. We did not find difference in our results using other statistical distributions 
of noise. We observe modulation instability with the fastest growing mode corresponding to 
wavelength 27r that develops to time-shifts t ^ 10 generating 16 peaks in such a box and leading 
to formation of one-dimensional wave turbulence. 

For our initial conditions kinetic energy is comparable to the potential one i^4, that means 
we are working in the regime of solitonic turbulence when solitons or quasi-solitons play significant 
role in the turbulent re-distribution of energy inside the system (see also [TT1|12]). In the integrable 
case Q the turbulence is called integrable and relaxes to one of infinite possible stationary states. 
Eq. (pl describes statistically irreversible movement to it's statistical attractor - one big soliton 
in the field of small perturbations that contain all the necessary information for time-reversal 
t ^ —t (compare to pTHTS] ). In case ofEq. (|3j) wave field ^ quickly (tens of nonlinear lengths) 
approaches to statistically steady state when high waves appear randomly in space and time. 

In our numerical simulations we used the 2nd-order Split-Step method in which linear and 
nonlinear parts of the equations were calculated separately. In order to improve simulations 
and save computational resources we employed adaptive change of spacial grid size Ax reducing 
it when Fourier components of solution at high wave numbers k exceeded 10~^^ max 
and increasing Ax when this criterion allowed. In order to prevent appearance of numerical 
instabilities, time step At also changed with Ax as At = hAx'^ with h < 0.1 (see [14^). 

Implementation of such numerical schema allowed us to safely shift up to t ~ 100 nonlinear 
lengths, i.e. when 

\\^hit,x)-^0.5h{t.x)\\/\\'^0.5hit.x)\\ < 1, (5) 

where ^/i(t,x) is the numerical solution calculated with a fixed coefficient h = At/Ax"^. We also 
did comparisons with 4th-order Split-Step method [T5[ [16] as well as with 4th- and 5th-order 
Runge-Kutta methods that confirmed validity of our results for single simulations up to t ~ 100 
nonlinear lengths. 

Thus, all methods we used, including comparisons with the results obtained with higher number 
of Fourier modes (lower Ax) with the same initial noise, gave virtually the same results for the 
classical NLS equation up to time shifts t ~ 60. Beyond 60 nonlinear lengths all methods and all 
schema parameters At and Ax we used gave different results in the sense of criterion ([5|. This 
behavior is connected with quasi-periodical dynamics of the classical integrable NLS equation: 
near t ^ 60 wave field again can be represented as ^ = l+^(x) with \(,{x)\ <C 1, but this time ^(x) 
contains numerical errors that are unique for numerical method and it's parameters At and Ax. 
Therefore, beyond 60 nonlinear lengths modulation instability develops differently for different 
numerical methods and their parameters. 

Nevertheless, comparison of our statistical results calculated with ho ~ 1/12, in particular 
PDFs, revealed no difference with the results obtained with h = 0.25ho or with the help of 
4th-order Split-Step or 4th- and 5th-order Runge-Kutta methods far beyond 100 nonlinear 
lengths. 

3. We start this section with our results concerning evolution of kinetic and potential energy for 
Eq. ([2]) -([3]) as well as for evolution of total energy and wave action for Eq. ([3|. As shown on FIG. 
la, generalized NLS equation accounting for six- and eight-wave interactions ([2| demonstrates 
asymptotic movement to its statistical attractor - one big soliton containing all the potential 
energy H4 + Hq + and immersed in the field of small fiuctuations. Up to time shifts t ~ 300 
for a = 0.04, P = 0.001 and t ~ 150 for a = 0.064, P = 0.002 six- Hq and eight- wave interactions 
energy are small compared to kinetic Hd and four- wave interactions energy H4. Beyond these 
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time shifts are seen the first signs of relaxation when Hd^ H^^ Hq and Hs become comparable 
first and then approach to their asymptotic values after t ~ 2000 for a = 0.04, (3 = 0.001 and 
t - 1500 for a = 0.064, (3 = 0.002. 



d) 
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FIG. 1: (Color on-line) Dependence on time shift t of averaged over ensemble (a) kinetic energy (Hd) 
(black), four- {Ha) (blue), six- {Hq) (green) and eight-wave interactions energy {Hg) (red) for generalized 
NLS equation accounting for six- and eight-wave interactions M), a = 0.040^ /3 = 0.001; (b) total energy 
{Hd -\- H4 -\- Hq), (c) wave action {N) and (d) kinetic energy{Hd) (black), four- {H4) (blue) and six- 
wave interactions energy {Hq) (green) for generalized NLS equation accounting for six-wave interactions, 
dumping and pumping terms Hb, a = 0.128^ di = 0.04^ dn = 0.0004^ p — 0.05. 

Another way of regularization of collapses due to six- wave interactions with the help of dumping 
and pumping terms in Eq. (|3| allows one to constantly stay in the regime where |i^6| <C |^4|, 
i.e. when dynamics of the system is close to that described by the classical NLS equation (see 
FIG. Ic). After the first 20-50 nonlinear lengths the system reaches statistically steady state when 
energy drain due to dissipation is compensated by energy income due to pumping and wave action, 
momentum and total energy as well as kinetic H^^ four- and six- wave interactions energy H/^ 
fluctuate near their mean values (see FIG. lb,c). The statistically steady state is determined 
by coefficients a, di^ dn and p only and does not depend on initial distribution ^|t=o: thus, we 
observe the same results for ^|t=o = e(x) where \e{x)\ <C 1 is stochastic noise. In the steady 
state spectra //e, spatial correlation functions g{x)^ correlation length Xcorr defined as full width 
at half maximum for g{x) and PDFs for wave field amplitude |^| measured at some specific time 
t no longer depend on time that allows us to perform additional averaging over time. 

FIG. 2 shows evolution of spectra, spacial correlation function and PDF for generalized NLS 
equation accounting for six- and eight- wave interactions ([2|. During relaxation process spectra 
gradually accepts some universal bell-shaped form with two exponential parts corresponding to 
small and large wave-numbers k and can be divided by two parts - the upper one corresponding 
to growing soliton and the lower one corresponding to small fiuctuations field. At small time- 
shifts t ~ 50 and for different sets of a and /3 normalized spatial correlation function g{x)/g{0) 
is close to universal form which approaches to Gaussian for small \x\ < Xcorr- At larger time- 
shifts universal form of correlation functions persists and becomes sharp, almost triangular in its 
maximum at x = 0. Evolution of correlation length Xcorr for Eq. ([2|-(|3| is shown on FIG. 3. 

Optical rogue waves are typically observed in the initial stage of modulation instability [TJ |5] 
that occurs at time shifts t ^ 10 for our initial pulse and noise amplitude. At these time 
shifts even despite significantly nonlinear regime ~ H4 when weak wave turbulence approach 
does not work PDFs for Eq. ([2| and for the classical NLS equation fiuctuate near Gaussian 
distribution. These fiuctuations occur because different realizations of noise lead to different 
times of modulation instability development and diminish at larger times (see FIG. 4a for classical 
NLS equation for example). 

Thus, and at t = 40 PDFs for Eq. ([2| turn out to be almost Gaussian with the exception of 
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small amplitudes regions. We observe such behavior until six- Hq and eight- wave interactions 
Hs are small compared with dispersion and four- wave interactions i^4, i.e. until dynamics 
of the system remains close to that described by the classical NLS equation. When Hd, i?4, 
Hq and become comparable there appears power-law region at middle amplitudes on the 
PDFs that divides initial exponential decay at low amplitudes and the peak at high amplitudes 
corresponding to relaxation to one big soliton. Beyond this peak PDFs again decay exponentially 
that in our opinion must be connected with deviations in amplitude of this big soliton due to 
deviations in initial data ^|^=o- Spectrograms enlarged at pulse maximums shown on FIG. 5a, 5b 
demonstrate that typical large wave events for Eq. ([2| are collisions of quasi-solitons at moderate 
time shifts t ~ 100 and single quasi-solitons at large time shifts t ~ 2000 where relaxation to 
statistical attractor occurs (see [8( IT7]). 

In contrast to Eq. Q-Q, Eq. (|3| demonstrates significantly different evolution of PDFs with 
time when long non- Gaussian tails appear already during modulation instability development 
t ~ 10. Beyond these time shits PDFs fiuctuate preserving their shape as shown on FIG. 
4b where PDFs are represented versus and is the averaged over ensemble 

squared amplitude at the corresponding time shift. At t ~ 50 fiuct nations of wave action and 
(I^P) diminish (compare with FIG. Ic) as the system reaches statistically steady state where 
energy drain due to dumping terms is compensated by the energy input from the pumping term. 
The corresponding "fat tail" in the region of large amplitudes decays faster than any power of 
amplitude |^| and turns out to be similar to const *exp(—|^|) as shown on FIG. 4c. Spectrogram 
on FIG. 5c demonstrates that a typical collapse in case ofEq. ^ originates as a collision of 
several quasi-solitons. 

We now compare results obtained for Eq. (|3| with different six-wave interactions coefficients a 
and dumping and pumping parameters di^ dn and p. First, we did several numerical simulations 
for a set of different di , dn and p with fixed a and found our results qualitatively do not depend 
on pumping and dumping parameters. Then we fixed di = 0.04, dn = 0.0004 and p = 0.05 
and did five simulations for different six-wave interactions coefficient from a = to a = 0.256. 
Corresponding spectra, spacial correlation functions and PDFs are shown on FIG. 6. Thus, 
in the statistically steady state spectra consists of two exponential parts ~ ex.-p{—q\k\) with 
increasing coefficient q from small to large wave-numbers k and a transition region between them 
corresponding to medium wave- numbers. Normalized spacial correlation functions g{x)/g{0) 
approach to some universal function of x/xcorr that for small x is close to Gaussian. 

The mean wave action at large time shifts changes with a from ~ 170 at a = to A/" ~ 70 
at a = 0.256 because the main power drain from the system occurs during regularization of 
collapses [9 , while the number of collapses per time unit depends on a: thus, six- wave interactions 
dominate for a = 0.032 starting from |^| > 6 and for a = 0.256 starting from |^| > 2. In the 
framework of Eq. ([3| there is no independent on a universal spatio-temporal dynamics for the 
collapses saturation. Since mean wave action at statistically steady state depends on six-wave 
interactions coefficient, mean square amplitude also depend on it: decreases with increasing 

a from 1.7 at a = to 0.7 at a = 0.256. 

Therefore on FIG. 6c we plot normalized PDFs depending on that allows us to 

examine PDFs for different a, di, dn and p all on one graph. At small amplitudes |^| PDFs 
demonstrate almost Gaussian decay that changes at moderate and high amplitudes to strongly 
non-Gaussian ~ exp(— 1^|) for a > 0, while for a = PDF remains Gaussian even for high 
waves. The latter in the sense of Gaussian PDF for the classical NLS equation ([T]) shows that 
dumping and pumping terms in Eq. (|3| do not affect the form of PDF. For non-zeroth a the 
corresponding non-Gaussian addition to PDF increase with a, i.e. frequency of large wave 
occurrence increase with six- wave interactions coefficient. 

4. We would like to underline two our main results. First, presence of nonlinearity and signifi- 
cantly nonlinear regime of a system do not necessarily mean non-Gaussian PDFs as demonstrated 
by our results for the classical NLS equation (IT]) shown on FIG. 4a and for generalized NLS equa- 
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tion accounting for pumping and dumping terms ([3]) with a = as shown on FIG. 6c. Generahzed 
NLS equation accounting for six- and eight-wave interactions ([2| also demonstrates Gaussian tails 
on the PDFs in the region of large amplitudes if its dynamics is close to that described by the 
classical NLS equation. 

Second, for Eq. ([2| we discovered emergence of power-law tails starting from moderate time- 
shifts when first signs of relaxation to statistical attractor appear. We also demonstrated presence 
of non-Gaussian tails ~ exp( — 1^|) when higher waves appear more frequently for generalized NLS 
equation accounting for six-wave interactions, pumping and dumping terms (|3| even when six- 
wave interactions are small compared to four- wave interactions. The corresponding non- Gaussian 
addition does not qualitatively depend on dumping and pumping parameters, disappear for zeroth 
six- wave interactions coefficient a = and increase with a. 

D. Agafontsev thanks E. Kuznetsov for valuable discussions concerning this publication, M. 
Fedoruk for access to and V. Kalyuzhny for assistance in author's work with Novosibirsk Su- 
percomputer Center. This work was done in the framework of Russian Federation Government 
Grant (contract no. 0035 with Ministry of Education and Science of RF, November 25, 2010), 
and also supported by the program of Presidium of RAS "Fundamental problems of nonlinear 
dynamics", program of support for leading scientific schools of Russian Federation and RFBR 
grant 09-01-00631-a. 



[1] D.R. Solli, C. Ropers, P. Koonath, B. Jalali, Optical rogue waves, Nature 450, 1054-1057 (2007). 
[2] J.M. Dudley, C. Finot, G. Millot, J. Gamier, G. Genty, D. Agafontsev and F. Bias, Extreme events 

in optics: Challenges of the MANUREVA project, Eur. Phys. J. Special Topics, Vol. 185 (2010), 

pp.125 - 133. 

[3] D.R. Solli, C. Ropers, and B. Jalali, Active Control of Rogue Waves for Stimulated Supercontinuum 
Generation, Phys. Rev. Lett., vol. 101, iss. 23 (2008), p. 233902. 

[4] J.M. Dudley, G. Genty, B.J. Eggleton, Harnessing and control of optical rogue waves in supercon- 
tinuum generation, Optics Express, Vol. 16, No. 6, pp. 3644-3651 (2008). 

[5] G. Genty, CM. de Sterke, O. Bang, F. Dias, N. Akhmediev, J.M. Dudley, Collisions and turbulence 
in optical rogue wave formation. Physics Letters A, vol. 374, iss 7 (2010), pp. 989 - 996. 

[6] A. Mussot, A. Kudlinski, M. Kolobov, E. Louvergneaux, M. Douay, and M. Taki, Observation of 
extreme temporal events in CW-pumped supercontinuum. Optics Express, vol. 17, iss. 19 (2009), 
p.17010. 

[7] K. Hammani, B. Kibler, C. Finot, A. Picozzi, Emergence of rogue waves from optical turbulence. 

Physics Letters A 374, 34 (2010) pp. 3585-3589. 
[8] D.S. Agafontsev, On the modulation instability development in optical fiber systems, 

arXiv:1004.3209vAl, Pis'ma v ZhETF, vol.91 (2010), iss. 12, pp. 699-704. 
[9] Y. Chung and P.M. Lushnikov, Strong Collapse Turbulence in Quintic Nonlinear Schrodinger Equa- 
tion, Physical Review E, v. 84, 036602 (2011). 
[10] W.H. Renninger, A. Chong, and F.W. Wise, Dissipative solitons in normal- dispersion fiber lasers, 

Phys. Rev. A77, 023814 (2008). 
[11] V.E. Zakharov, A.N. Pushkarev, V.F. Shvetz, V.V. Yan'kov, Solitonic turbulence, Pis'ma v Zh. Eksp. 

Teor. Fiz. 48, 79-81 (1988) [JETP Lett. 48, 83-85 (1988)]. 
[12] V. Zakharov, F.Dias, A. Pushkarev, One- Dimensional Wave Turbulence, Physics Reports, 398, 1, 
1-65, (2004). 

[13] R.Jordan, C.Josserand, Self- organization in nonlinear wave turbulence, Phys. Rev. E 61, 1527 - 1539 
(2000). 

[14] T.L Lakoba, Stability analysis of the split-step Fourier method on thebackground of a soliton of the 
nonlinear Schrodinger equation, arXiv: 1008.4974vl 

[15] G.M. Muslu, H.A. Erbay, Higher-order split-step Fourier schemes for the generalized nonlinear 
Schrodinger equation. Mathematics and Computers in Simulation, v. 67, iss. 6 (2005). 

[16] R.L Mclachlan, On the numerical integration of ordinary differential equations by symmetric com- 
position methods, SIAM Journal on Scientific Computing, v. 16, iss. 1 (1995). 



7 



[17] E.B. Treacy, Measurement and Interpretation of Dynamic Spectrograms of Picosecond Light Pulses, 
J. Appl. Phys.42, 3848-3858 (1971). 



8 



T = 40 



T = 400 



T = 2000 







1 








-100 c 
h 

T = 


) 100 
40 










, , >\ 



-4 -2 2 4 

x/x 

corr 



T=40 



10^ 

10-^ 

10-^' 
1 

0.8 
0.6 
0.4 
0.2 






a) 











-100 
k 

T = 400 



100 









b) 











4-2024 

x/x 

corr 



T=400 



10^ 

10-^ 

10-^ 
1 

0.8 
0.6 
0.4 
0.2 








A 


a): 




i 







-100 




k 

T = 2000 



100 













, rb) 










/' 

: y 1 


. .»y 



4-2024 

x/x 

corr 



T=2000 




20 40 
T=40 



20 40 
T=400 



20 40 
T=2000 




10^ 



10" 



10" 



10' 




FIG. 2: (Color on-line) Solid lines are spectra Ik (a), normalized spatial correlation functions g{x)/g{0) 
(h) and PDFs in semi-log (c) and log-log scale (d) for generalized NLS equation accounting for six- and 
eight-wave interactions ^ for a — 0.040^ (3 — 0.001 (red) and a — 0.064^ /3 = 0.002 (blue) at time 
shifts t — 4Q, t — 400 and t — 2000. Dashed lines for graphs (h) and (c) are Gaussian distributions, for 
graphs (d) - power-law tails. Linear parts on graphs (d) demonstrate power-law decay and continue from 
= 3 = 20. 
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FIG. 3: Correlation length Xcorr defined as full width at half maximum for spatial correlation function g{x) 
depending on time t for (a) generalized NLS equation accounting for six- and eight-wave interactions 
a = 0.040^ /3 = 0.001 and (h) generalized NLS equation accounting for six-wave interactions, dumping 
and pumping terms a = 0.256, di = 0.04, dn ^ 0.0004, p = 0.05. 
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FIG. 4: (Color on-line) PDFs for (a) classical NSL equation at time shifts t = 10 (red), t — 12 (blue), 
t = 30 (black) versus \^\^ , (b) generalized NLS equation accounting for six-wave interactions, dumping 
and pumping terms Q at time shifts t — Q (purple), t = 8 (red), t = 10 (blue), t — 12 (greenjand 
in the statistically steady state averaged over time [50, 250] (black) versus |^|^/(|^|^), (c) Eq. \^ in 
the statistically steady state averaged over time [50, 250] versus |^|. Solid lines - PDFs, dashed lines - 
Caussian distributions for graphs (a) and (b) anc? exp(— 1^|) tail for graph (c). Craphs (b) and (c) are 
calculated with parameters di — 0.04, dn — 0.0004, p — 0.05, a — 0.128. 
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FIG. 5: ( Color on-line) Spectrograms of typical large wave events for generalized NLS equation accounting 
for six- and eight-wave interactions gj, a = 0.040, P = 0.001, at moderate (a) \^\max = 6.1, t = 72.7 
and large (b) time shifts \^\max — 7.7, t — 1789 and (c) generalized NLS equation accounting for six-wave 
interactions, dumping and pumping terms \^\max — 7.9 at t — 219.7. 
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FIG. 6: (Color on-line) Solid lines are (a) spectra Ik, (h) normalized spatial correlation functions 
g{x)/g{0) and (c) normalized PDFs for generalized NLS equation accounting for six-wave interactions, 
dumping and pumping terms ^ with fixed di = 0.04, dn = 0.0004, p = 0.05 and a = (purple), 
a — 0.032 (red), a = 0.064 (blue), a = 0.128 (green) and a = 0.256 (black) in the statistically steady 
states. Dashed lines for graphs (b) and (c) - Gaussian distributions. 



